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ABSTRACT 

The Zeeman pattern of Mn I lines is sensitive to hyperfine structure (HPS) and, because of this reason, they 
respond to hG magnetic field strengths differently from the lines commonly used in solar magnetometry. This 
peculiarity has been employed to measure magnetic field strengths in quiet Sun regions. However, the methods 
applied so far assume the magnetic field to be constant in the resolution element. The assumption is clearly 
insufficient to describe the complex quiet Sun magnetic fields, biasing the results of the measurements. The 
diagnostic potential of Mn I lines can be fully exploited only after understanding the sense and the magnitude 
of such bias. We present the first syntheses of Mn I lines in realistic quiet Sun model atmospheres. Plasmas 
varying in magnetic field strength, magnetic field direction, and velocity, contribute to the synthetic polarization 
signals. The syntheses show how the Mn I lines weaken with increasing field strength. In particular, kG 
magnetic concentrations produce Mn I A5538 circular polarization signals (Stokes V) which can be up to two 
orders of magnitude smaller than what the weak magnetic field approximation predicts. As corollaries of this 
result, (1) the polarization emerging from an atmosphere having weak and strong fields is biased towards the 
weak fields, and (2) HPS features characteristic of weak fields show up even when the magnetic flux and 
energy are dominated by kG fields. For the HFS feature of Mn I A5538 to disappear the filling factor of kG 
fields has to be larger than the filling factor of sub-kG fields. Since the Mn I lines are usually weak. Stokes V 
depends on magnetic field inclination according to the simple consine law, a scaling that holds independently 
of the magnetic field strength. Atmospheres with unresolved velocities produce very asymmetric line profiles, 
which cannot be reproduced by simple one-component model atmospheres. Inversion techniques incorporating 
complex magnetic atmospheres must be implemented for a proper diagnosis. Using the HFS constants available 
in the literature we reproduce the observed line profiles of nine lines with varied HFS patterns. Consequently, 
the uncertainty of the HFS constants do not seem to limit the use of Mn I lines for magnetometry. 

Subject headings: line: profiles - Sun: magnetic fields - Sun: photosphere 



1. INTRODUCTION 

When the polarimetric sensitivity and the angular resolu- 
tion exceed the required threshold, magnetic signals show 
up almost everywhere on the solar photosphere. The sig- 
nals in supergranulation cell interiors are particularly weak, 
however, most of the solar surface is in the form of cell in- 
teriors and, therefore, these weak signals probably set the 
total unsigned magnetic flux and magnetic energy exist- 
ing in the photosphere at any given time (e.g., jUnnoi 1959 ; 
Stenflol [T982L lYi et all Il99l ISanchez Almeida! [19981 12004 : 
Schrijver & Titld l2003h . The importance of these ubiq 



iTrujillo Bueno et aTl |2004|) . Measuring implies assuming 
many properties of the unresolved complex magnetic field 
and, in doing so, the measurements become model depen- 
dent and prone to bias. Depending on the technique used 
for m easuring, the real quiet Sun exh ibits weak daG fields 
(e^ lStenfloril982l: lFaurobert-Scholll il993: Bommieret aD 



20051). intermediate hG fields (e.g., ' Lin & Rimmele. .19991: 
Khomenko et al. | [ 2003 : lEopez Ariste e t all 120061) . or strong 



uitous quiet Sun magnetic fields depends, to a large ex- 
tent, on the magnetic field strengths characterizing them. 
For example, the magnetic flux and the magnetic energy 
scale as powers of the field strength, and the connec- 
tivity between photosphere and corona is probably pro- 
vided by the magnetic concentrations with the largest field 
strengths (e.g., Dommguez Cerdena et al. 2006a). Unfor- 
tunately, measuring quiet Sun magnetic field strengths is 
not a trivial task. All measurements are based on the 
residual polarization left when a magnetic field of com- 
plex topology is observ ed with finite angular resolution 
(e.g.. lEmonet & Cattanea2 001: S anchez Almeida et al.l2003l: 
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kG fields (e.g., 'Sanchez Almeida & LitesI (200tf: 'LitesI I200; 
Ipommguez Cerderia et al. 2003). Such discrepancies can 
be naturally understood if the true quiet Sun contains a 
continuous distribution of field strengths going all the way 
from zero to two kG. Different techniques are biased differ- 
ently and, therefore, they tend to pick out a particular part 
of the distribution. This scenario is very much consistent 
with realistic numerical simulations of magneto-convection 
fcatta neo 1999; Stein & Nordlund 2006 : Vogler et al. 200"1 
IVogle r & Schiissler 2007j). In order to provide a compre- 
hensive observational description of the quiet Sun magnetic 
field strengths, one has to combine different methods care- 
fully chosen to have complementary biases. Then the full 
distribution can be assambled tak ing the biases into account 
(iDominguez Cerdena et al.ll2006al) . 

In an effort to complement t he existing magnet i c field 
^tteitgthniiMg»Q*nB@£6fiaaifipaK, iLopez Ariste et all (|2002|) 
proposed the use of spectral lines whose Zeeman patterns 
are sensitive to hyperfine structure (HFS). The formal- 
ism to deal with the HFS of spectral lines in magnetic 
atmospheres was developed more than thirty years ago by 
iLandi Degl'Innocentil (Il975h . According to such formalism, 
the polarization of the HFS lines vary with magnetic field 
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strength very differently from the hnes commonly used in 
sol ar magnetometry. This u nusual behavior was invoked 
by iLopez Ariste et all (|2002|) when proposing the use of 
HFS Mn I lines as a diagnostic tool for magnetic field 
strengths. Lopez Ariste and coworkers have applied the idea 
to measure magne tic field stre ngths in quiet Sun region s 
(ILopez Ariste et al.l l2006. 200j; lAsensio Ramos et al.ll2Bo7h . 
Since the number of observables is limited, they minimize 
the statistical error of the measuremet by minimizing the 
number of free parameters to be tuned. Milne-Eddington 
atmospheres (ME) are used to synthesize the polarization of 
the lines. The magnetic field is assumed to be constant and, 
therefore, the measurements provide some kind of weighted 
average of the true field strengths existing in the resolution 
elements. As we pointed out above, the topology of the quiet 
sun magnetic field is complex, with a distribution of field 
strengths and polarities coexisting in a typical resolution 
element. Then the ill-defined average provided by the Mn I 
lines is expected to be biased toward a particular range 
of field strengths, as it happens with the rest of the quiet 
Sun measurements (e.g., the NIR Fe I lines overlook kG 
magnetic c oncentrations, Sanchez Almeida & Lited l2000l 
ISocas-Navarro & Sanch ez Almeida 2003; the traditional 
visible Fe I lines exagge r ate the c ontribution of kG fields, 
Bellot Rubio & CoUadosI I2003L iMartmez Gonzalez et all 
2001 " the Hanle scatteri ng polarization signals are not sen- 
sitive to hG and kG, Faurobert et al.l l200ll ISanchez Almeidal 
12005). The bias presented by the HFS Mn I lines is so far 
unknown, and the existing and forthcoming measurements 
based on those lines will be fully appreciated only when the 
sources of systematic effects are properly acknowledged and 
quantified. In order to explore the magnitude and the sense of 
the expected effects, we undertake the synthesis of Mn I lines 
in a number of realistic quiet Sun atmospheres with complex 
magnetic field distributions. The main trends are presented 
here. 

The paper is organized as follows. Section |2] describes the 
software developed to carry out the syntheses. First, a ME 
code is needed to compare our syntheses with the results in the 
literature. Then, a plane parallel one-dimensional (ID) code 
allows us to explore the influence of realistic thermodynamic 
conditions on the polarization of lines with HFS. Finally, a 
Micro- Structured Magnetic Atmosphere (MISMA) code pro- 
vides additional realism to the modeling since it includes cou- 
pling between magnetic field strengths and thermodynamic 
conditions, magnetic fields varying along and across the line- 
of-sight (LOS), mixed polarities in the resolution element, etc. 
All these c odes are based on the origina l FORTRAN routine 
written bv lLandi Degrinnocentil (119781) . The analysis is fo- 
cused on the line most often used in observations, namely, 
MnlA5538. We discuss its intensity and circular polariza- 
tion, since the linear polarization signals are very weak and 
they remain undetected so far Single component and multi- 
component syntheses of this line are described in §|3]and §|4l 
respectively. An exploratory attempt to consider unresolved 
mixed polarities and velocity fields is carried out in § [5] The 
polarization of other Mn I lines is also considered in § [6l The 
implications of these syntheses are discussed in § |7] 



2. DESCRIPTION OF THE SYNTHESIS PROCEDURES 

The synthesis of Stokes profiles has been carried out using 
three different approaches to solve the radiative transfer equa- 



tion for polarized light, 

dl 

---K(I-S), (1) 

as 

where I = {I,Q,U,V)^ is the column vector containing the 
Stokes parameters, S = (Bt, 0,0,0)^ is the source vector with 
Bj the LTE source function, K is the 4x4 absorption ma- 
trix and, finally, the variable of integration s corresponds to 
the length along the LOS. ( For details on matrix elemen ts and 
the rest of the equation, see'Landi Degrinn ocentill976i ) The 
three codes differ in the model atmosphere used to describe 
the photosphere. The first one solves equation ([T]) under ME 
hypothesis. The source function is assumed to vary linearly 
with the continuum optical depth r. Then the radiative trans- 
fer equation admits an analytic solution involving the two co- 
efficients of the source function {Bt -Bq+Bit) and the seven 
coefficients of the absorption matrix, wh ich are assumed to be 
constant (see lLandi Degrinnocentill975 : equation [14]). The 
second code solves equation ([TJ for one-dimensional (ID) at- 
mospheres. Realistic temperatures, pressures and magnetic 
field vectors are considered, but they only vary along the LOS 
so that a single ray suffices to synthesize the spectrum. The 
third code solves a different version of equation ([T]). It is 
obtai ned by a spatial averag ing under the MISMA hypoth- 
esis ( Sanchez Almeida et al.|[l996) . It assumes the solar at- 
mosphere to have inhomogeneities smaller than the photon 
mean-free -path whose details are washed out by the radiative 
transfer process. Only their average properties matter. The 
new equation to be solved is, 

f =-(K)(I-S'), (2) 

as 

where S' = (K) ' (KS) is a mean source vector. The brack- 
ets in equation ^ denote average over a volume whose 
dimension along the LOS is of the order of the photon mean- 
free -path. As a result of the averaging, (K) and S' are slowly 
varying functions of s, with this slow variation accounting 
for the large scale structure in the atmosphere. Pursuing 
the idea that optically-thin details are irrelevant for the final 
spectrum, Sanchez Almeidal (11997) puts forward a type of 
model MISMA made of distinct components with different 
temperatures, pressures and magnetic field vectors. The coex- 
istence of different components imposes a number of physical 
constraints to be satisfied, in particular, the lateral pressure 
balance causes the magnetic field to be coupled with the 
thermodynamics (see the original reference for an exhaustive 
description). These model MISMAs, with their constraints, 
provide the degree of realism required to carry out our ex- 
ploratory work. Note that the model MISMAs reproduce all 
kinds of Stokes profiles observed in the quiet Su n, including 
those with extreme asymmet ries (Sanchez Alm eida & LitesI 
120001; 'Socas-Navarro & San chez Almeida l2002t 
Dommguez Cerderia et al. 2006b). Equation ([T]) is for- 
mally identical to equation (|2]). They are integrated using a 
predictor-corrector method with the variables equi-spaced in 
a fixed grid of atmospheric heights. As we mention in the 
introduction, the linear polarization signals Q and U are very 
weak, and we do not analyze them here despite the fact that 
they are synthesized together with / and V. 

All the codes are based on the FORTRAN procedure HY- 
PER by Landi Degl'Innocenti (1978) that computes the Zee- 
man pattern of any line with hyperfine structure. The pat- 
tern depends on the externally imposed magnetic field, to- 
gether with a set of atomic parameters, namely, the hyper- 
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fine structure constants, the quantum numbers of the two lev- 
els involved in the transition, the relative isotopic abundance, 
and the isotope shifts. When the Zeeman splitting is com- 
parable to the hyperfine structure splitting, the Stokes profiles 
are strongly disrupted showing typical hyperfine structure fea- 
tures. The pattern depends on the so-called hyperfine structure 
constants A and B, which account for the two first terms of 
the Hamiltonian describing the interaction between the elec- 
trons in an atomic level and the nuclear magnetic moment - 
A corresponds to the magnetic-dipole coupling whereas B de- 
scribes the electric-quadrupole coupling. The constants A and 
B can be determined empirically or theoretically. The values 
of A adopted for our Mn I syntheses are listed in Table [T] The 
table includes the appropriate references, which indicate that 
B is negligible small in all cases. 

2. 1 . Testing the codes 

We carry out different types of tests using the line 
Mn I A5538 as target. Intensity profiles are compared with the 
unpolarized solar atlas. The syntheses are also compared with 
the (ME) Stokes profiles existing in the literature. The codes 
are cross-checked one with respect to another, and with exist- 
ing synthesis codes that neglect the HPS. All these efforts are 
summarized next. 

The first series of tests is based on non-magnetic atmo- 
sp heres. We start by reco vering the ME Stokes / profile 
in iLopez Ariste et"aD (l2002l Fig. 3) when using the appropri- 
ate ME parameters (Lopez Ariste 2007, private communica- 
tion). In addition. Fi g. [T] sho ws a reasonable ME fit to the 
FTS solar spectrum (iNeckell 1199 9) using a Doppler width 
AAd = 0.029 mA, a damping constant a - 0.3, an absorp- 
tion coefficient 770 = 2.7, a source function given by Bq = 1 -0 
and B\ = 1.1, and the atomic parameters in Table [T] Figure [T] 
also includes a fit to the observed profile using the ID code. 
The line is sy nthesized in the quiet Sun model atmosphere by 
iMaltbv et al.l (IT986h using the atomic parameters in Table [T] 
The fit was carried out by trial-and-error varying the velocity 
along the LOS, the macroturbulence (v,„ac, set to 1.3 km s~'), 
and the Mn abundance (log[Mn], set to 5.30 in the scale where 
log[H]=12). The MISMA code was also tested against the 
observed Stokes / assuming the thermodynamic parameters 
of the model atm osphere to be given by the model quiet Sun 
in ISanchez Alm eida (1997, § 4.1). As in the case of ID at- 
mospheres, the atomic parameters are those listed in Table [T] 
and we vary the velocity along the LOS, the macroturbulence 
{vmac =1-5 km s~'), and the Mn abundance (log[Mn]=5.31). 
One additional consistency test was carried out with the ID 
code. We compare the synthesis obtained neglecting the hy- 
perfine structure (i.e., imposing A - B - 0) with the intensity 
of the ID code from which we inherited the absorption coef- 
ficient routines (iSanchez Almeidall 19921 § 5). The profiles of 
the two syntheses are identical within the numerical precision. 

The second series of tests refers to magnetic atmospheres. 
Employing the atomic parameters and the atmospheres de- 
scribed above, we include a constant magnetic field vector. 
We use three different magnetic field strengths, B - 100 G, 
300 G, and 900 G, with the same inclination, 9-60°, and 
azimuth, (f) = 0°. The ME syntheses reproduce the Stokes pro- 
files in Fig. 4 of lL6pez"A riste et al. (2002). Stokes V profiles 
synthesized with the three codes are shown in Fig. [2l The 
agreement between the syntheses supports their internal con- 
sistency - if the codes are tuned to produce similar Stokes / 
profiles (Fig.fTJ, then for the same field strength they also pro- 
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Fig. 1. — Mm A5538 Stokes / profiles syntliesized using tiie three codes 
and compared with the profile from the solar FTS atlas. The inset identifies 
the different codes and the observation. The wavelengths, in units of A, are 
referred to the laboratory wavelength of the line. The profiles are normalized 
to the quiet Sun continuum intensity. 




Fig. 2. — Mm A5538 Stokes V profiles synthesized using the three codes 
with B = 100 G (a) and B = 900 G (b). The insets pair each type of line with 
the corresponding synthesis code. The profiles are normalized to the quiet 
Sun continuum intensity, and the wavelengths are refereed to the laboratory 
wavelength of the line. 

duce similar Stokes V (Fig. |2ll. The Stokes V profiles repre- 
sented in Fig. [2] correspond to (and are c onsistent with) the 
two behaviors of Mn I A5538 described by Lopez Aris te et al.l 
([2002). The left panel in Fig. [2] shows the kind of profile 
for B < 700 G, with a characteristic reversal not far from the 
line core (hereafter HFS hump). Figure [Jh corresponds to 
B > 700 G, which do not generate HFS humps at the core. The 
presence or absence of such reversal is used by Lopez Ariste 
and co-workers to distinguish between sub-kG fields (present) 
and kG fields (absent). 

The results in § [6] can be regarded as a third indepen- 
dent test. The same non-magnetic quiet Sun model MISMA 
that reproduces the observed intensity of MnlA5538 also 
accounts for eight additional Mn I lines with very different 
HFS patterns. Such agreement shows the realism and self- 
consistency of the syntheses. 

3. SYNTHESES OF Mn I A5538 IN ID ATMOSPHERES 

The weak magnetic field approximation is routinely used 
in solar magnetometry to, e.g., calibrate magnetograms. It 
is also valid for lines with HFS (iLandi Deal' In nocentill 19751: 
Landi Degl'Innocenti & Landolfi 2004), and so it was used 
by Lopez Ariste et al. (2006) to estimate the relative contri- 
bution of weak and strong fields to the observed quiet Sun 
signals. When the approximation holds the Stokes V profile 
is proportional to the longitudinal component of the magnetic 
field BcosO, 

V = -c ^Bcos9, (3) 
dX 

where the constant c depends on the particu- 
lar spectral line, and the derivative of the inten- 
sity profile, dl/d\, gives the variation with wave- 
length of Stokes dLandiJegrinnocentil Il975t 
iLandi Degl'Innocenti & Landolfil |2004|) . We tested the 

Note, incidentally, that given the observed intensity profile of 
Mm A5538 (Fig. [T}, a hump at its Stokes V line core is unavoidable when 
B^O. 
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TABLE 1 

Atomic parameters used in the synthesis of Mn i lines 



Line 


Wavelength^ 


Transition^ 




^Oggf glow" 










[A] 




[eV] 






[10"^cm-i] 


[10^3 cm-i] 



Mn I A4457 
Mn I A4470 
Mn I A4502 
Mn I A4739 
Mn I A5420 
MniA5516 
MniA5538 
MnlA8741 
Mm A 15262 



4457.04 
4470.14 
4502.22 
4739.11 
5420.36 
5516.77 
5537.76 
8740.93 
15262.70S 



6p0 



■3/2 



^ "3/2 



'3/2" 
''5/2" 



«'03/2- 

/2 



y ^5/2 
y ^3/2 

,.6p0 

^3/2 



9/2 



7/2 



5/2_ 



3.073 
2.940 
2.919 
2.940 
2.142 
2.178 
2.186 
4.434 
4.889 



-0.555" 
-0.444" 
-0.344" 
-0.490" 
-1.460" 
-1.847" 
-1.920"= 
-0.05 f 
+0.45' 



1.875 
1.198 
1.368 
1.198 
1.584 
1.866 
3.327 
1.712 
2.000 



1.759 
1.200 
1.427 
0.400 
1.886 
2.400 
2.400 
1.554 
2.2848 



15.61' 
1.7<= 
-4.6<= 
1.7<= 
13.5'' 
16.5'' 
27.8'' 
-13'' 
25.2*' 



22.8'' 
6.4'= 
1.5'= 
22.3<= 
-21.5'' 
-31.5'' 
-31.5'' 
15.5'' 
27.5'' 



Note. — The symbols undefined in the main text stand for: lower-level oscillator strength gf, lower-level Lande 
factor giow> upper-level Lande factor gup, and lower-level excitation potential x- 

" NIST Atomic Spectra Database, 'Ralche nko eT^ (20071) '' ILefebvre et 5] <2003l) '= IBlackwell- Whitehead etal] 
120051) '' Ifisher&Peck (1939)"= Margrave. <1972D '' Set to match the observed solar spectrum with 
log[Mn]=5.3 l.g .Asensio Ramos et al.i t2007i) 
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Fig. 3. — Maximum Stokes V amplitude as a function of the magnetic 
field strength for ME, ID, and MISMA model atmospheres. The prediction 
corresponding to the weak field approximation is represented as the dotted 
line. All magnetic field inclinations are set to 60° . Magnetic field strengths 
are given in G. 

approximation with the quiet Sun model atmospheres de- 
scribed in § 12.11 which are realistic in the sense that they 
reproduce the observed intensity profile when the magnetic 
field strength is close to zero (Fig. [T]). The syntheses assume 
the magnetic field to be constant along the LOS, with the 
inclination 6 set to 60°. Figure [3] shows the variation of the 
maximum Stokes V signal as a function of the magnetic field 
strength. The weak field approximation breaks down for 
fairly low magnetic field strengths, i.e., B > 400 G (compare 
the dashed line and the dotted line in Fig. |3]l. Even more, 
the deviations are very important for kG magnetic field 
strengths. When the field strength is 1.5 kG the weak field 
approximation yields a Stokes V signal twice larger than the 
synthetic signals in ID model atmospheres (Fig. [5]l. In order 
to test the dependence on cos 6* predicted by equation ([3]), 
we also carried out syntheses with various magnetic field 
inclinations. This dependence turns out to be closely followed 
by the synthetic signals, even for strong kG magnetic field 
strengths. The behavior was to be expected since Mn I A5538 
is a weak line satisfying the weakly polarizing medium 
approximation (Sanchez Almeida & Truiillo Bueno 1999). In 
this approximation Stokes V scales with its specific absorp- 
tion coefficient, which is proportional to cos 6 independently 
of whether or no the spectral line has HFS. 

The quiet Sun contains plasmas with all magnetic field 
strengths from to 2 kG (§ [TJ. According to a mechanism 
originally put forward by ISpruili (1 1976ft and now well estab- 
lished, one expects a strong correlation between the magnetic 
field strength and the temperature of the observable photo- 
spheric layers. Plasmas having kG fields must be strongly 




Fig. 4. — Stokes / and V profiles synthesized in three ID model atmo- 
spheres corresponding to the quiet Sun (the solid line), a network region (the 
dotted dashed fine), and a plage region (the dashed line). The magnetic field 
is constant with a strength of 100 G and an inclination of 60°. AH profiles 
have been normalized to the quiet Sun continuum intensity. The dotted line 
is just a reference to show the continuum level. 

evacuated to stay in mechanical balance within the photo- 
sphere. Consequently, strongly magnetized plasmas are trans- 
parent, showing light coming from deep and therefore hot 
(sub-)photospheric layers. In order to explore the dependence 
of the Mm A5538 Stokes V signals on the thermodynamics, 
we compare synthetic profiles formed under different ther- 
modynamic conditions but the same magnetic field strength. 
The line is synthesized in a quiet Sun model atmo sphere 
jMaltbvetal. 1986), a network model atmosphere ( Solankil 
[1986; Solanki & Steenbock.l988j. and a plage mo del atmo- 
sphere (Solankil 1 19861: ISolanki & Steenbocklll988h . with the 
temperature increasing from quiet Sun to network. The re- 
sults are shown in Fig.|4] First note that the continuum inten- 
sity of Stokes / increases from quiet Sun (the coolest model) 
to network (the hotest model). In anticorrelation with contin- 
uum intensity, the Stokes V signals decrease with increasing 
temperature (Fig. |4l right). This behavior is not attributable 
to the HFS, but it is due to the ionization balance of Mn in 
the photosphere. Mn I is the minor species so that a small 
modification of the ionization equilibrium does not alter the 
Mn II abundance but it drastically changes the number of Mn I 
atoms available for absorption. Increasing the temperature in- 
creases the ionization, and the Mn I lines weaken. Figure [5] 
shows the relative abundance of Mn I as a function of height 
in the atmosphere for a constant density representative of the 
layers where the lines are formed (say, 100 km above con- 
tinuum optical depth equals one). According to Fig. [S] if the 
presence of a strong magnetic field decreases the formation 
height of a spectral line by 100 km, then the Mn I is depleted 
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Fig. 5. — Relative abundance of Mn I and Mn 11 as a function of the height 
in the atmosphere for a quiet Sun stratification of temperatures. (The temper- 
ature corresponding to each height appears on the top of the figure.) The gas 
pressure is the same for all the heights, and it has been chosen to match the gas 
pressure of the non-magnetic quiet Sun some 100 km above the continuum 
optical depth equals one. The Fe 1 abundance is included for comparison - its 
overabundance with respect to Mn I is due to its higher ionization potential. 

by more than one order of magnitude (compare the values at 
100 km and km). Note that this behavior is common to all 
Mn I lines, and it is also typical of other minor species like 
Fe I (e.g., the dotted line in Fig.jS]). In short, the Mn I lines are 
expected to weaken with increasing magnetic field strength. 

The coupling between magnetic field strength and temper- 
ature is fully accounted for in the MISMA syntheses, which 
partly explains the dimming of the kG Stokes V signals shown 
in Fig.|3](the solid line). Another part of the dimming is due to 
non-magnetic plasma obscuring the magnetic plasma, and ef- 
fect sp ecific to MISMAs and elaborated by .Sanchez Almeidal 
(l2000h . 

4. MULTICOMPONENT SYNTHESIS OF Mn I A5538 

The syntheses carried out so far assume very simple mag- 
netic configurations. The magnetic field is either constant 
or, in the case of model MISMAs, one magnetic component 
varies along the LOS. Due to the complex topology of the 
quiet Sun magnetic fields, our limited angular resolution, and 
the spatial averages carried out to obtain appropriate signals, 
those simplifications are insufficient to represent observed 
Stokes profiles. Real polarization signals are formed in plas- 
mas having a range of magnetic properties, therefore, we im- 
prove the realism of the synthesis assuming multi-component 
model atmospheres. Specifically, we assume 

/= / '"'^ P(B)IsdB, (4) 
Jo 

and ^ 

I {co^e)P{B)VBdB. (5) 

The symbols Ib and Vb represent the Stokes / and V profiles 
produced by an atmosphere with a single longitudinal mag- 
netic field of strength B, whereas P{B) stands for the fraction 
of resolution element occupied by such atmosphere. P{B) is 
usually referred to as magnetic field strength probability den- 
sity function (PDF). It is normalized to one, and it can be 
envisaged as the filling factor corresponding to each B. The 
maximum field strength existing in the quiet Sun sets B^ax 
to some 2 kG. As we show in Appendix [A] equations Q 
and ((Sj hold even for an arbitrary distribution of magnetic 
field directions, provided that the average inclination factor 
(cos 9) is computed properly. We would like to make it clear, 
however, that the PDF approach for representing the compli- 
cations of the quiet Sun magnetic fields also has limitations. 
Equations Q and (O imply a one-to-one relationship between 
the magnetic field and the thermodynamic structure of the at- 
mosphere. Although a coupling between magnetic field and 
thermodynamics is to be expected, the real relationship should 
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Fig. 6. — Magnetic field strength PDFs used for multi-component synthe- 
ses. The one labelled DC is from Dommsuez Cerdeiia et al. 1 2006a) while 
the VIO, V50 and V200 are from iVoglen t2003l) . See text for details. 
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Fig. 7. — Multi-component Stokes V profiles synthesized using ME atmo- 
spheres (a), ID atmospheres (b), and MISMAs (c). The labels and the lines 
in the insets refer to the PDFs in Fig. [6] The profiles are normalized to the 
quiet Sun continuum intensity. 

have a significant scatter (e.g. IVogleij|2003h . Magnetic field 
variations along the LOS are restricted, and the overlapping 
along the LOS of structures with different magnetic fields is 
ignored. Our PDF approach is only a first approximation to 
the problem. 

The four PDFs in Fig. |6] are used for synthesis. They cor- 
re spond to one of the semi-ernpirical quiet Sun PDF obtained 
by lDomfnguez Cerdena et all (l2006al) (labelled DC), plus the 
P DFs of the three magneto-convection numerical simulations 
bv^Vogled (l2003h . representing magnetic fluxes of 10 G, 50 G 
and 200 G (VIO, V50 and V200, respectively). When ME 
atmospheres, ID atmospheres, and model MISMAs are com- 
bined according to these PDFs, one obtains the Stokes V pro- 
files shown in Fig. |7] In this case the magnetic field inclina- 
tion is constant and set to 60°, i.e., (cos0) = 1 /2. The corre- 
sponding Stokes / profiles for the ME and the ID cases are 
practically identical to those in Fig. [T] The MISMA synthe- 
sis has a continuum intensity 5% lower than that of the quiet 
Sun, since the thermodynamic structure of the model MISMA 
used for synthesis is cooler than the typical quiet Sun model 
atmospheres (see Fig.l 1 in Sanchez Almeida & Lites 2000). 

Note that the Stokes V profiles in Figs.|7]always show HFS 
humps. This result warns us against simplistic interpretations 
of the presence of HFS features as the unequivocal signature 
of an atmosphere dominated by hG magnetic fields. Despite 
the fact that most of the atmospheric volume has very low 
field strengths, the magnetic flux and the magnetic energy of 
some of the PDFs used for synthesis are dominated by the tail 
of kG fields. The extreme case corresponds to V200 in Fig. [6l 
where 60% of the magnetic flux and 87% of the energy^ is in 

^ Fo llowing Sanchez Almeidal <2004t) and Dommguez Cerd ena et alj 
(2006i), the unsigned magnetic flux and the magnetic energy are computed 
as the two first moments of the PDF. 



6 



Fig. 8. — Relative contribution to the average Stokes V profile of a multi- 
component atmosphere characterized by the PDF labelled as DC in Fig.|6] 
The integral of the function represented in the plot is proportional to the 
Stokes V signal at a particular wavelength. The three curves correspond to 
any wavelength under the weak field approximation (the solid line), the wave- 
length where the average Stokes V signals is maximum (the dashed line), and 
the wavelengths where the HFS hump shows up (the dotted line). The units 
of the ordinates are relative since the curves have been scaled to the solid line 
when B — > 0. 

field strengths larger than 700 G, which is the field strength 
where the HFS hump disappears according to the ME syn- 
theses. Even in this case the Stokes V profiles present clear 
HFS humps (see Fig. |71 the triple dotted-dashed lines). In 
order to understand this behavior one has to realize that the 
various magnetic field strengths existing in the atmosphere do 
not contribute to the observable signal in proportion to their 
magnetic flux, as expected if the weak magnetic field approxi- 
mation would be satisfied (equation lO). The strong kG fields 
are always under-represented (§ [3]), and this bias is particu- 
larly severe at line center where the HFS hump shows up. The 
effect is shown in Fig. [8] which contains the integrand of equa- 
tion ^ for one of the PDFs in Fig. [6] (DC). The integrand is 
given in three cases, (a) assuming the weak field approxima- 
tion to hold (the solid line), (b) at the wavelength where the 
average Stokes V is maximum (the dashed line), and (c) at 
the line core (the dotted line). The individual spectra have 
been synthesized in model MISMAs, but the qualitative be- 
havior is common to all other model atmospheres. According 
to Fig.[8]the contribution of the kG fields to the HFS hump of 
the average Stokes V signal is some ten times smaller than that 
predicted by the weak field approximation. The relative con- 
tribution of the various field strengths to the hump is worked 
out as part of § 14.11 

4. 1 . On the diagnostic provided by the HFS hump 

The observations of MnlA5538 analyzed in the literature 
employ single component ME atmospheres (see §[T]). The ob- 
served Stokes V profiles are fitted with synthetic profiles to 
obtain mean magnetic field strengths. In order to explore the 
diagnostic content of such inversion technique when applied 
to complex magnetic atmospheres, we developed a simple 
least squares minimization procedure to fit single-component 
ME profiles Vb to the synthetic multi-component V. The pro- 
cedure fits 

V = aBVB, (6) 

where the factor ub and the field strength B are the only two 
unknowns. Figure |9] shows the single-component ME fit of 
the multi-component Stokes V produced by the DC PDF. In 
the case that the syntheses are based on ME atmospheres 
(Fig. |9t), the match is excellent. If the synthesis is based 
on model MISMAs, the fit worsens a bit but still the agree- 
ment is fair (Fig.|9]5). In both cases the field strength resulting 
from the ME inversion is about 550 G. We have unsuccess- 
fully tried to link this average magnetic field strength with 
a particular property of the underlying PDF (the solid line 
in Fig. (6]). It is neither the unsigned flux (~150 G) nor the 
most probable field strength (~13 G). Figure[8l the solid line. 




Fig. 9. — Least squares fits of Stokes V profiles synthesized in multi- 
component atmospheres (the solid lines) using single-component ME atmo- 
spheres (the symbols), (a) The multi-component Stokes V profile is synthe- 
sized in ME atmospheres, (b) The multi-component Stokes V profile is syn- 
thesized in model MISMAs. The values for ag and B retrieved from the fits 
are included in the figures. The two synthetic profiles correspond to the same 
PDF (DC). The profiles are normalized to the quiet Sun continuum intensity. 

shows BP(B), i.e., the unsigned flux per unit of magnetic field 
strength. It peaks at 50 G and 1700 G which, again, do not 
seem to bear any obvious relationship with the average ME 
magnetic field strength (~550 G). 

In order to understand the relationship between the HFS 
hump and the magnetic fields existing in the multi-component 
atmosphere, we adopt a strategy avoiding ME fits. The pres- 
ence or absence of a hump is certainly associated with the 
balance between weak fields and strong fields, but it is un- 
clear how the hump is related with the fraction of atmosphere 
filled by weak and strong fields. One can quantify the pres- 
ence of the hump using the derivative of the Stokes V profile 
at a wavelength slightly bluewards of the HFS maximum 
(i.e., centered in linear growth that ends up in the hump; see 
Fig|2^ where c^O). Then the sign of 



dV_ 



(7) 



indicates the presence (> 0) or absence (< 0) of a hump^. 
According to equation (O, 



{cos 9) P{B) 



dVB 
dX 



dB, 



(8) 



so the slope of the mean Stokes V is just a weighted mean of 
the slopes that characterize each field strength. Consequently, 
dVB/dX\x^. yields the contribution of each field strength to the 
sign of the mean profile or, in other words, it yields the con- 
tribution of each field strength to the presence or absence of a 
hump in the average profile. 

We have applied equation dH) to understand the HFS humps 
of the synthetic profiles in Fig. |71 Figure [TOb shows the 
Stokes V derivative at line core (Ac = 0) for the ME, the ID and 
the MISMA model atmospheres. The magnetic field strength 
where this derivative is zero indicates when Vg changes from 
having to laking of a hump. As the figure shows, this tran- 
sition field strength depends on the model atmosphere. The 
positive and negative lobes of the curves in Fig.fTOb have sim- 
ilar amplitudes, which indicates that for the average Stokes 
profile to lack of a hump, the filling factor of the kG fields in 
the atmosphere has to be similar to the filling factor in weak 
fields. The requirement is not fulfilled by the PDFs in Fig. [6l 
whose filling factors are dominated by weak fields. Such su- 
periority explains why all profiles in Fig.|7]show HFS humps. 
Equation (HI) also explains why the magnetic flux in the form 
of kG structures has to exceed by far the flux in weak fields for 
the HFS hump to disappear. The HFS hump goes away when 

^ This criterion assumes positive magnetic polarity, implying a positive 
Stokes V blue lobe. The same rule applies to negative polarity if one reverses 
the sense of the inequalities. 



Fig. 10. — (a) Sensitivity of the HFS hump to the distribution of magnetic 
field strengths in the resolution element. The hump appears or goes away 
depending on the integral of dVB/d\\x^ weighted with the PDF. The three 
curves correspond to the three model atmospheres used for multi-component 
synthesis (see the inset), (b) The derivative shown in (a) weighted with the 
V200 PDF. The integral is positive meaning that the resulting Stokes V pro- 
files have HFS humps; see the triple dotted-dashed lines in Fig. [7] 

weak fields and strong fields have similar filling factors, but 
then the kG fields dominate the flux since the (unsigned) mag- 
netic flux scales as the filling factor times the field strength. 
For example, consider a resolution element filled half-and- 
half with structures of 100 G and 1 kG. The HFS hump is 
not present, but the magnetic flux in kG fields is ten times the 
weak field flux. 

Among the three curves in Fig.fTOb. the one corresponding 
to MISMAs reflects the solar conditions best, since it includes 
the dimming of the line with increasing field strength. This 
curve has positive and negative lobes of approximately the 
same area, therefore, the presence-absence of HFS humps in 
the Mn I A5538 Stokes V profiles is controlled by the fraction 
of photospheric atmosphere with field strengths smaller-larger 
than some 1 kG (i.e., the zero crossing point of the curve). In- 
terpreting real observations with this rule-of-thumb requires 
disregarding the atmospheric volume having very weak mag- 
netic fields, which produce no polarization and therefore do 
not contribute to the observed signal (see FigfTOb). 

5. MIXED POLARITIES AND UNRESOLVED 
VELOCITIES 

We expect the quiet Sun resolution elements t o con- 
tain mixed polarities (ISanchez Almeida & LitesI l2000t 
ISocas-Navarro & Sanchez Almeidal 120021) . which leads to a 
reduction of the Stokes V signals. In the multi-component 
approximation parlance (equation JSj), having unresolved 
mixed polarities implies, 

(cos6') < (|cos6'|), (9) 

since positive and negative cos 6 tend to cancel when carry- 
ing out the averages. If the positive polarity {cosO > 0) and 
the negative polarity (cos < 0) have the same distribution of 
inclinations but they fill different fractions of resolution ele- 
ment, then 

(cos0) = (|cos0|)(2/-l), (10) 

with / the fraction filled by the magnetic fields of posi- 
tive polarity. If / does not depend on the magnetic field 
strength, the effect of mixed polarities is only a global scal- 
ing of the Stokes V profile with no effect on Stokes I (equa- 
tions m, ISl, and ifTOl ). However, the above view of the 
effect produced by unresolved mixed polarities is too sim- 
plistic. It implicitly neglects the existence of unresolved 
velocities in the resolution element coupled with changes 
of polarity. When the sense of the velocity and the polar- 
ity are correlated, the Stokes V of opposite polarities are 
Doppler shifted. They do not perfectly cancel, leaving ob- 
servable residuals. Such coupling has been observed in Fe 
lines, and it produces very asymmetric Stokes V profiles 



Fig. 1 1 . — Examples of Stokes V profiles synthesized in semi-empirical 
model MISMAs having unresolved mixed polarities correlated with veloci- 
ties. They correspond to different classes in Sanchez Al meida & Litesi >200Q) 
- (a) class 6, (b) class 9, and (c) class 7. (d) Offspring of the class 7 model 
atmosphere where the two polarities have the same sign and the velocity dif- 
ference has been increased. The profiles are normalized to the quiet Sun 
continuum intensity. 




Fig. 12. — Two component Stokes V synthesis of weak and strong fields 
having different Doppler shifts. The strong fields are blueshifted with respect 
to the weak fields - see (a). The sum of the profiles in (a) leads to the profiles 
in (b), the solid line. If the relative Doppler shift is not considered, then 
the average Stokes V corresponds to the dotted line in (b). The profiles are 
normalized to the quiet Sun continuum intensity. 

(ISanchez Almeida & LiteslbOO tf). In order to illustrate the ef- 
fect. Fig. [TT] shows Stokes V profiles of MnlA5538 synthe- 
sized in empirical model MISMAs with mix ed polarities and 
different velocities (ISanchez Almeida & Litesii20 00). 

Velocities and other magnetic quantities are also ex- 
pected to be coupled. For example, strong kG fields ap- 
pear in intergranular lanes with downflows , whereas weak 
fields prefer granules with upflows (e.g., ICattaneol Il999t 
ISocas-Navarro et 120041) . This combination of Doppler 
shifts and magnetic field strengths renders profiles like the two 
component average shown in Fig.[T2b. the solid line. It corre- 
sponds to the sum of the profiles in Fig.[T2h. which are shifted 
by some 3.5 km s~^ The HFS hump now appears elevated 
from the zero level, and the Stokes V blue wing is broadened 
with a conspicuous double aborption structure. These two 
features are commo n among the o bserved Stokes V profiles 
of Mm A5538 (Lope z Ariste et al.l f2006. Figs. 8 and 9). The 
same features appear in the MISMA synthesis of Fig.fTTB. 

In short, mixed polarities and velocity gradients are prop- 
erties of the real quiet Sun magnetic fields that affect the ob- 
served polarization in a significant way. They are not easy to 
deal with but must be taken into account for a realistic inter- 
pretation of the observations. 

6. SYNTHESIS OF OTHER Mn I LINES 

The purpose of synthesizing Mm lines different from 
Mn I A5538 is threefold. First, it confirms that our LTE syn- 
theses reproduce all kinds of observed HFS patterns. Second, 
it offers a chance to explore whether other Mn I lines show a 
sensitivity to magnetic fields complementary to Mm A5538. 
Finally, it proves that the uncertainty of the HFS constants do 
not limit the use of Mn I lines for magnetometry. 

All the syntheses in this section employ the non-magnetic 
quiet Sun MISMA that reproduces the Mn I A5538 profile in 
Fig. [T] Varying the macroturbulence and, in some cases also 
tuning the oscillator strength, we managed to fit the observed 
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Fig. 13. — Mm lines without obvious HFS features in the observed solar 
spectrum. The solid lines show our syntheses whereas the dotted lines corre- 
spond to the solar atlas lNeckel'1999). Wavelengths, given in A, are refereed 
to the laboratory wavelength of the lines. Each line is identified by the wave- 
length in the label of the ordinate axis. The profiles are normalized to the 
quiet Sun continuum intensity. 




Fig. 14. — Mm lines with HFS features. The solid lines correspond to 
our syntheses whereas the dotted lines show the solar atlas (Neckel 1999). 
Wavelengths, given in A, are refereed to the laboratory wavelength of the 
lines. Each line is identified by the wavelength in the label of the ordinate 
axis. The profiles are normalized to the quiet Sun continuum intensity. 

Stokes / profiles of all the lines we tried^. A total of eight 
Mn I lines reproduced by the code are shown in Figs. [T3] and 
[T4l They are split in two sets; Fig. [13] includes lines without 
obvious HFS features in the observed solar spectrum, whereas 
Fig.[T4]shows lines with varied and assorted HFS. The agree- 
ment between the synthetic profiles and the solar atlas is good, 
an impression reinforced by the fact that the HFS constants 
determining the HFS patterns have not been tuned, but they 
come directly from the literature (Table[T]). 

The dependence of the various lines on magnetic field 
strengths has been studied by synthesizing them in the mag- 
netic quiet Sun model MISMA used in §|4l increasing the field 
strength from ~ 100 G to ~ 2000 G. All the lines share a com- 
mon behavior - they weaken with increasing field strength, 
and the HFS features tend to disappear. This behavior is also 
sh ared by the near IR line Mm A 15262 recently pointed out 
bv lAsensio Ramos et all (120071) - Fig. bottom right. On 
top of this general trend, and depending on the actual HFS 
constants, the Stokes V profiles of different lines present spe- 
cific features. The line Mn I A55 16 behaves like Mn I A5538, 
with a single HFS reversal at the line core. Then it is possible 
to relate the presence or absence of the hump with the dis- 
tribution of field strengths in the atmosphere, exactly as the 
analysis carried out in 34.11 The patterns of Mn I A5420 and 
Mn I A8741 are more complicated, and we have not been able 
to agree on a particular feature of the Stokes V profiles that 
can be treated as in 94.11 

7. DISCUSION AND CONCLUSIONS 

The Zeeman pattern of the Mn I lines depends on hyper- 
fine structure (HFS), which confers them a sensitivity to hG 
magnetic field strengths different from the lines traditionally 

' The pattern of Mn I A5407 puzzled us for some time. As suggested by R. 
Manso, it turned out to be a blend of the Mn I line with two fairly strong Fe I 
lines that remain unidentified in the classical list by Moore et al. ( 1966), and 
which do not appear in the NIST Atomic Spectra Database. We found them 
with approximately the observed strengths among the list of lines by Kurucz 
at |http : //kurucz . h arvard .edu/linelists . htmlj 



used in solar magnetometry. This peculiarity has been used 
to measure magnetic field strengths in quiet Sun regions (see 
§ [T]). The methods applied so far assume the magnetic field 
strength to be constant in the resolution element, an approxi- 
mation driven by feasibility rather than based on physical or 
observational arguments. Actually, it is not a good approxi- 
mation since the magnetic fields of the quiet Sun are expected 
to vary on very small spatial scales, with field strengths span- 
ning from zero to 2kG. Under these extreme conditions, all di- 
agnostic techniques employed in quiet Sun magnetometry are 
strongly biased towards a particular range of field strengths, 
and the Mn I signals are not expected to be the exception. 
Consequently, the diagnostic content of the Mn I lines cannot 
be fully exploited unless their biases are properly understood. 
Such task is undertaken in the paper by exploring the response 
of Mn I lines in a number of realistic quiet Sun scenarios. 

Three complementary LTE synthesis codes have been writ- 
ten and tested (ME, ID and MISMA; § They provide a 
number of relevant results, the first one being the ability to 
reproduce all observed Mn I HFS patterns. We reproduce the 
observed unpolarized line profile of nine assorted lines corre- 
sponding to all HFS sensitivities (§[6]). The study is focused 
on the line most often used in magnetometry, MnlA5538, 
however its behavior should be representative of the other 
lines. According to the weak magnetic field approximation, 
the Stokes V signals scale with the longitudinal component 
of the magnetic field, i.e., the magnetic field strength times 
the cosine of the magnetic field inclination with respect to 
the LOS. We verify that the scaling on cosine holds very 
tightly. The scaling on the magnetic field strength, how- 
ever, breaks down soon (when B > 400 G for Mm A5538). 
Even for ME atmospheres, the weak field approximation pre- 
dicts a Stokes V signal twice as large as the synthetic one for 
B ~1.5 kG. When the expected coupling between the thermo- 
dynamic conditions and the magnetic field strengths is taken 
into account, the dimming of the kG Stokes V signals can be 
as large as two orders of magnitude (see Figs. [3] and [Sj. The 
dimming of the polarization signals formed in kG magnetic 
concentrations affects all Mn I lines (§ [3]), and it has signifi- 
cant observational implications. 

If the resolution elements contains both weak and strong 
fields, then the kG fields tend to be under-represented in 
the average profile. We have modelled the bias assuming a 
multi component atmosphere, where the synthetic signals are 
weighted means of the Stokes profiles corresponding to each 
single field strength. The weight is given by the fraction of 
atmosphere filled by each field strength, i.e., by the magnetic 
field strength probability density function PDF (equations lH 
and iQ). According to our modeling, even when the (un- 
signed) magnetic flux and the magnetic energy are dominated 
by kG magnetic fields, the Stokes V profile of MnlA5538 
can show HFS reversal at line core characteristic of hG fields 
(Fig- El). A pure morphological inspection of the Stokes pro- 
files does not suffice to infer which is the dominant magnetic 
field strength in the resolution element. For the HFS hump 
of MnlA5538 to disappear the kG filling factor has to be 
larger than the sub-kG filling factor and, consequently, when 
the HFS hump disappears the magnetic flux and magnetic en- 
ergy of the atmosphere are completely dominated by kG fields 

(§133]). 

Detecting Stokes V profiles with HFS features indicates 
the presence of hG fields in the resolution element. How- 
ever, this sole observation does not tell whether the hG 
field strengths dominate. There seem to be two extreme al- 
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tematives to exploit the diagnostic potential of these lines. 
First, improving the spatial resolution of the observations to 
a point where the quiet Sun magnetic structures can be re- 
garded as spatially resolved. Unfortunately, this possibil- 
ity does not seem to feasible at present. Realistic simu- 
lations of magneto-convection indicate that quiet Sun mag- 
netic fields are uniform only at the diffusive l ength scales 
(e.g. ICattaneol \l99^ IVogler & Schussleii 120071) . which are 
of the order of a few km in the photosphere (e.g. ISchiissleJ 
119861) . These scales are very far from the angular resolution 
of the present measurements, and even much smaller than 
the le ngth-scale for the radiative t ransfer average along the 
LOS (ISanchez Almeida et al.lll996l) . We prefer the alternative 
possibility, namely, developing inversion techniques where 
complex magnetic atmospheres are included into the diagnos- 



tics. Using appropriate constraints, the number of free pa- 
rameters required to describe such atmospheres can be main- 
tained within reasonabl e limits (e.g., the model MISMAs in 
ISanchez Almeidalll997l) . Dealing with unresolved velocities 
also favors detailled inversion codes (§[5]). We are presently 
working on these improvemen ts needed to devel op the diag- 
nostic technique pioneered by lLopez Ariste et all 
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APPENDIX 

MULTI-COMPONENT ATMOSPHERE WITH MAGNETIC FIELD VECTORS VARYING IN DIRECTION 

The equations ^ and ^ used for multi-component syntheses hold even when the magnetic field vector presents a distribution 
of magnetic field inclinations, 9, and magnetic field azimuths, 0. In the general case the PDF depends on B, 9 and 4>, V(B, 9, (f)), 
so that the average Stokes / and Stokes V profiles are. 



2n 



ViB, 9, (/)) S[iB, 9, (t)) dBd9d(j), 



(Al) 



2-IT 



ViB, 9, 4>)Sv{B, 9, 4>) dBd9d(l), 



(A2) 



with V(B, 9, (f)) normalized to one. The symbols Si(B, 9, 0) and Sv(B, 9, 4>) stand for the Stokes / and V profiles corresponding to 
each magnetic field vector. As the tests carried out in §[3] show, and as it is to be expected from the weakly polarizing medium 



10 



approximation (ISanchez Almeida & T rujillo Bu enol|l999l) . the intensity of the weak Mn I lines is independent of inclination and 
azimuth so that Sj(B, 6,(f)) c:^ Ig. Similarly, Stokes V is independent of and depends on 9 through the factor cos 9, Sv(B, 6,cf))- 
VbCos9. The symbol Vg represents the Stokes V profile for longitudinal magnetic field (cos6' = 1). Then from equations (lAll ) and 
(IA2I ) one retrieves the equations Q and ^ used in the main text, 



/= / P{B)hdB, (A3) 
Jo 

V= I {cos 9) PiBWedB, (A4) 



^0 

with 

r27V rTV 

P{B)^ / / V{B,9,cb)d9d(j), (A5) 

and 

(cos6') = / / cos9ViB,9,(t))d9d(t). (A6) 

PiB) Jo Jo 



